Estimating hitting probabilities of an interacting particle system on a graph

نویسنده

  • T. Brereton
چکیده

We consider the problem of estimating hitting probabilities related to a class of interacting particle systems. These systems, in which two types of particles — ‘electrons’ and ‘holes’ — move on a graph, are simplified versions of models describing charge transport in disordered materials. The probability of interest is the probability that an electron reaches a certain region of the graph before colliding with a hole. We provide a detailed description of our model and explain how it can be simulated. Next, we give a brief introduction to importance sampling, which we use to improve the efficiency of our estimators. To our knowledge, importance sampling has not yet been used to estimate probabilities related to interacting particle systems. We describe how importance sampling can be used to improve the efficiency of estimators of hitting time probabilities involving discrete time Markov chains. We then use importance sampling to estimate the probability that we are interested in. In doing so, we observe that there are a number of complexities that arise when working with interacting particle systems. We describe some simple heuristics for implementing importance sampling. These heuristics make minimal changes to the probability measure under which the original system evolves. We consider a specific example of our problem and investigate the effectiveness of the importance sampling approach. We show that our estimators outperform standard Monte Carlo estimators. Finally, we describe possible future work, which includes a more sophisticated importance sampling approach that uses ‘locally optimal’ changes of measure.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

A functional central limit theorem for interacting particle systems on transitive graphs

A nite range interacting particle system on a transitive graph is considered. Assuming that the dynamics and the initial measure are invariant, the normalized empirical distribution process converges in distribution to a centered di usion process. As an application, a central limit theorem for certain hitting times, interpreted as failure times of a coherent system in reliability, is derived.

متن کامل

Genetic genealogical models in rare event analysis

We present in this article a genetic type interacting particle systems algorithm and a genealogical model for estimating a class of rare events arising in physics and network analysis. We represent the distribution of a Markov process hitting a rare target in terms of a Feynman–Kac model in path space. We show how these branching particle models described in previous works can be used to estima...

متن کامل

Sequential Monte Carlo simulation for the estimation of small reachability probabilities for stochastic hybrid systems

The problem of estimating the probability that a system reaches a given set within some time horizon is considered. Standard Monte Carlo methods for reachability probability estimation do not require specific assumptions on the system under consideration. However, they are computationally demanding when the probability to be estimated is small. An Interacting Particle System (IPS) approach has ...

متن کامل

Rare Event Simulation Using Reversible Shaking Transformations

We introduce random transformations called reversible shaking transformations which we use to design two schemes for estimating rare event probability. One is based on interacting particle systems (IPS) and the other on time-average on a single path (POP) using ergodic theorem. We discuss their convergence rates and provide numerical experiments including continuous stochastic processes and jum...

متن کامل

Mean-Field Limit of a Stochastic Particle System Smoothly Interacting Through Threshold Hitting-Times and Applications to Neural Networks with Dendritic Component

In this article we study the convergence of a stochastic particle system that interacts through threshold hitting times towards a novel equation of McKeanVlasov type. The particle system is motivated by an original model for the behavior of a network of neurons, in which a classical noisy integrate-and-fire model is coupled with a cable equation to describe the dendritic structure of each neuron.

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2015